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= (54) Title: METHOD FOR BIOMECHANICAL SIMULATION OF A SET OF OSSEOUS JOINTS 

== (54) Titre : PROCEDE DE SIMULATION BIOMECANIQUE DUN ENSEMBLE D 'ARTICULATIONS OSSEUSES 

^5 (57) Abstract: The invention relates to a method for biomechanical simulation of a set of osseous joints of a patient, in particular 
jj^^ rachis. The inventive method consists in recording a digital three dimensional model embodied at least partially in the form of rigid 
= bodies interconnected by means of joints in a reference position, personalising the model geometry by specific data of the patient 
in said reference position, personalising said digital model by particulating interaction parameters [mobility or rigidity character- 
== istics] of each joint connecting said rigid bodies according to detected client's characteristics. Said method is characterised in that 
^= the particularisation of the interaction parameters consists in obtaining the space position of at least the part of the rigid bodies and 
555 interpolating for determining the calculated position of other rigid bodies in order to produce a numerical index containing the rela- 
tive position of each rigid body, performing at least one defined constraint on the patient and collecting information on the general 
balance position of said patient [other than the reference position], and in determining analytical functions which make it possible to 
approximate the interaction parameters, thereby reproducing the measured relative positions for each couple of rigid bodies. 
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(57) Abrcge : La prSsente invention concerne un precede de simulation biomecanique d'un ensemble d'articulations osseuses d'un 
patient, notamment du rachis, comportant une etape d'enregistrement d'un modele numerique tridimensionnel representee au moins 
en partie par des corps rigides relies par des articulations, dans une position de reference; une etape de personnalisation de la geo- 
metric du modele par des donnees specifiques a un patient dans ladite position de reference; une 6tape de personnalisation dudit 
modele numerique par particularisation des parametres d'interaction [de mobilites ou caracteristiques de rigidit£] de chacune des 
articulations reliant lesdits corps rigides en fonction des caracteristiques constatees sur le patient ; caracterise en ce que l'£tape de 
particularisation des parametres d'interaction consiste a acquerir les positions dans l'espace d'une partie au moins des corps rigides, 
et a proc6der a une interpolation pour determiner la position calculee des autres corps rigides pour construire une table numerique 
comportant les positions relatives de chacun des corps rigides ; exercer au moins une contrainte d6terminee sur le patient et a acquerir 
une information sur la position gen6rale d'equilibre r6sultante du patient [autre que la position de reference]; determiner des fonctions 
analytiques permettant d'approximer les parametres d'interaction afin de reproduire les positions relatives mesmrees, ce pour chaque 
couple de corps rigide. 
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PROCEDE DE SIMULATION B I OMECANI QUE D / UN ENSEMBLE 
D' ARTICULATIONS OSSEUSES 



La presente invention se rapporte au domaine des 
5 logiciels de simulation biomedicale. 

La presente invention se rapporte plus 
particulierement a un procede de simulation biomecanique 
d'un ensemble d' art iculat ions osseuses d'un patient, et 
notamment du rachis. Ce proc£d£ permet de calculer, estimer 
10 et visualiser les consequences d'une operation chirurgicale 
sur une articulation. Le systeme a ete mis au point pour 
les operations de la colonne vertebrale et en particulier 
les techniques de stabilisation de la colonne vertebrale. 
D'une maniere generale, le systeme permet d' informer le 
15 chirurgien sur l'etat de l'^quilibre et la repartition des 
efforts dans la colonne du patient avant et aprds la 
simulation. 



L'art ant^rieur connait deja, par la demande de 

20 brevet am§ricain US 5 995 738 , un dispositif et une methode 
pour faciliter 1 ' implantation de composants artificiels 
dans les articulations. L # invention decrit des dispositif s 
et des m^thodes destines a determiner une position 
d' implant pour au moins un composant artificiel dans une 

25 articulation et a faciliter son implantation. L' invention 
comprend la creation d'un modele de 1 ' articulation du 
patient et la creation d'un module du composant a 
implanter. Les modeles crees sont utilises pour simuler le 
mouvement de 1 ' articulation du patient selon la position du 

30 composant. Ce document de l'art anterieur propose done un 
moyen de simulation physique et non virtuel pour analyser 
les mouvements d'une articulation et d'un implant. 

II est propose, dans la demande de brevet americain 
US 6 205 411, un planif icateur de chirurgie assist§e par 

35 ordinateur et un systeme de guidage intra-opSrat if . 



WO 2004/061721 



2 



PCT/FR2003/003943 



L' invention concerne un appareil destine a faciliter 
r implantation d'un composant artificiel dans une 
articulation. L 7 appareil comprend un predicteur geometrique 
et un simulateur biomecanique du mouvement pr£-operatif , 
5 c'est-a-dire qu'une sSrie de simulations est effectuee sur 
1' implant et 1 ' art iculat ion avant de proceder a 
1' operation. 

II est propose, dans la demande de brevet americain 
10 US 5 62 5 577, une methode informatique d' analyse de 
mouvement utilisant la dynamique . L' invention decrit une 
methode permettant d' analyser et d'afficher les mouvements 
d'un etre humain. Le corps du sujet est divis£ en une 
pluralite de segments relies entre eux par des 
15 articulations. Une fois le corps modelise ainsi, il est 
possible de simuler et d' analyser le mouvement du sujet. 

II est egalement propose, dans la demande de brevet 
PCT WO 99/06960, un systeme et procede permettant de 

20 d6finir et d f utiliser des comportements destines £ des 
chaines articulees dans des animations par ordinateur. En 
1 1 occurrence, on definit pour la chaine articulee au moins 
une commande telle qu'une commande de forme ou une commande 
de plan de resolution, et on utilise cette commande comme 

25 une contrainte reprise par le moteur d' animation pour 
l 1 animation de la chaine avec une cinematique inverse. 
Chaque commande comporte au moins deux cles, chaque cle 
comprenant un couple constitu6 d'un vecteur de sens 
effecteur et d'une contrainte associ^e. Dans le cas des 

30 cles de commande de forme, les contraintes associees 
comprennent une liste d 1 orientations pr6f erentielles des 
membres du corps. Dans le cas des cles de commande de plan 
de resolution, les contraintes associees comportent une 
orientation pref erentielle du plan de resolution. Quels que 

35 soient les buts assignes, les cles de commande 



WO 2004/061721 



3 



PCT/FR2003/003943 



selectionn^es subissent une interpolation mettant en oeuvre 
des ponderations appropri^es visant £ l'obtention d'une 
contrainte resultante a faire utiliser par le moteur 
d 1 animation. 

5 L'art anterieur connait egalement, par la demande 

de brevet PCT WO 99/60939, un systeme chirurgical 
interactif assiste par ordinateur pour aider le chirurgien 
a positionner des implants dans le femur ou des vis dans 
les pEdicules des vertebres. Ces systemes apportent une 

10 assistance & la navigation. Cependant, c'est une assistance 
axee principalement sur le positionnement , le reperage et 
le guidage des ancillaires. Cette technique de navigation 
chirurgicale est assez r^pandue, mais ne permet pas de 
repondre a toutes les interrogations du chirurgien. En 

15 effet, 1 ' utilisation de ces systemes permet de d^finir la 
trajectoire optimale de la vis p^diculaire, mais ne permet 
pas de savoir par exemple, si elle a ete plac6e sur la 
bonne vertebre . 

Dans la publication scientifique « Morphometric 

20 evaluations of personalised 3D reconstructions and 
geometric models of the human spine » Aubin C-E et al - 
Medical and biological engineering and computing, Peter 
Perenigrus Ltd. Stevenage, GB - 01/11/1997, plusieurs 
methodes de reconstructions geom£triques sont evalu^es. 

25 Dans la conclusion de cette publication scientifique, les 
auteurs indiquent qu'une des methodes evaluees est 
optimale. L' etude du Dr Aubin confirme qu'une 
reconstruction geometrique d'une colonne vertebrale a 
partir de plusieurs radiographics constitue une approche 

30 pertinente et reproductible pour 1' Etude et 1' Evaluation 
des dysf onctionnements du rachis ainsi que pour sa 
modelisation par la methode des elements finis. Cette 
publication ne porte aucunement sur les aspects de 
personnalisation mecanique in vivo (mobility et 

35 anthropometrie) , ni sur la simulation de strategies 
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op^ratoires, ni sur leurs methodes d' opt imisation par 
approche biomgcanique . Enfin l'approche mise en aeuvre dans 
la presente invention associant plusieurs techniques 
(personnalisation geometrique, mobilite et anthropometric 
5 in vivo, analyse numerique, critdres d'equilibre, methodes 
de recalage entre donnees internes et externes, calculs des 
efforts) n'y est ni abordee ni suggeree. 

Enfin, l'art anterieur connait, par la publication 
scientifique « A MR I based semi-automatic modeling system 

10 for computational biomechanics simulation » Hayasaka T et 
al - Medical imaging and augmented reality, 2001, 
International workshop on 10-12 june 2001 - 10/06/2001, un 
precede de generation semi -automati que d'un maillage. Le 
procede d^crit dans ce document s' applique au systeme 

15 cardio-vasculaire, e'est-a-dire aux tissus mous et non aux 
articulations osseuses comme la presente invention. A 
partir d'un modele issu d'une base de donnees et d'une 
image IRM, un algorithme cr€e un modele sp^cifique en vue 
d'une etude (nous presumons de la tenue m^canique du tissu) 

20 base sur la methode des elements finis. L'auteur cite par 
ailleurs des diff icult6s restant a resoudre en vue de cr^er 
un maillage utilisable. La simulation m^canique des tissus 
mous qui reste a ce jour du domaine de la recherche 
f ondamentale, est confrontee k la difficulty de les 

25 caracteriser mecaniquement in vitro et encore plus in vivo. 
Cette etude ne porte pas sur la simulation d'un ensemble 
d' articulations osseuses. Elle porte sur un proced6 de 
generation semi-automatique d'un maillage. 

30 Ces deux publications scientif iques traitent de 

probl^matiques differentes de celle de la presente 
invention. En particulier, elles ne mentionnent nullement 
d'etape de personnalisation d'un modele numerique par 
particularisation des parametres d' interaction [de 

35 mobilites ou caracteristiques de rigidites] de chacune des 
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articulations reliant lesdits corps rigides en fonction des 
caract^ristiques constatees sur le patient. Cette etape est 
f ondamentale dans la presente invention. 

5 La presente invention entend remidier aux 

inconvenients de l'art anterieur en permettant de simuler 
une operation de correction locale ou globale de courbures 
de la colonne vert^brale ou de pose d' implant vertebral a 
partir d' images radiographiques , de series d' acquisitions 

10 de caract£ristiques mesurees in vivo sur le patient et 
d'une base de donn^es d' implants. La presente invention 
permet de plus de simuler l'etat de la colonne juste apres 
1' operation chirurgicale . 

Bien entendu, la presente invention ne s' applique 

15 pas uniquement au rachis. Elle s' applique ggalement a 
d'autres articulations osseuses comme le genou. 

L' invention consiste en un systeme chirurgical 
assiste par ordinateur, permettant au chirurgien de 

20 simuler, en preoperatoire, les effets sur le patient de la 
chirurgie correctrice qu'il envisage de mettre en 
application. Ce systeme lui permettant de simuler ainsi 
plusieurs strategies operatoires, il off re au chirurgien un 
outil d'aide au choix de la strategie operatoire off rant le 

25 meilleur compromis entre la stabilisation et la mobilite. 



A cet effet, la presente invention concerne dans 
son accept ion la plus g€n£rale un proc£d£ de simulation 
biomecanique d'un ensemble d' art iculations osseuses d'un 
30 patient, notamment du rachis, comportant : 

• une etape d' enregistrement d'un modele 
num6rique tridimensionnel representee au moins en 
partie par des corps rigides relies par des 
articulations, dans une position de reference ; 
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• une etape de personnalisat ion de la 
geometrie dudit modele [position relative dans 
1 ' espace de chacun desdits corps rigides] par des 
donnees specifiques a un patient dans ladite 

5 position de reference [par exemple des 

radiographies du patient] ; 

• une etape de personnalisation dudit modele 
num§rique par part icularisat ion des parametres 
d' interaction [de mobilites ou caracteristiques de 

10 rigidity] de chacune des articulations reliant 

lesdits corps rigides en fonction des 
caracteristiques constat6es sur le patient ; 



caracterise en ce que 

15 

l'§tape de part icularisat ion des parametres 
d' interaction [de mobilites ou caracteristiques de 
rigidity] consiste & : 



20 * acquerir les positions dans 1' espace d'une partie 

au moins des corps rigides, et & proceder & une 
interpolation pour determiner la position calculee 
des autres corps rigides pour construire une table 
num^rique comportant les positions relatives de 
25 chacun des corps rigides ; 

• exercer au moins une contrainte determinee 
sur le patient et a acquerir une information sur la 
position ggnerale d'£quilibre r^sultante du patient 
[autre que la position de reference] ; 
30 • determiner des fonctions analytiques 

permettant d'approximer les parametres 
d' interaction [de mobility's ou caracteristiques de 
rigidite] afin de reproduire les positions 
relatives mesurees [fonctions desdites contraintes 
35 appliquees sur le patient et des modifications 
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geometriques observ§es sur le patient, et resultant 
de ces contraintes] et ce pour chaque couple de 
corps rigide . 

5 De preference, le modele numerique est defini par 

des parametres de positions geometriques des corps rigides 
et par des parametres de rigidite des articulations reliant 
les corps rigides. 

Avantageusement , l'etape de representation du 

10 resultat d'une contrainte consiste a recalculer le modele 
personnalise [en position d'equilibre] resultant d'un 
ensemble de contraintes [par exemple implantation d'une 
prothese ou d'un implant] comprenant au moins une 
contrainte statique exercee sur au moins deux corps 

15 rigides, et imposant un posit ionnement relatif avec une 
mobilite ou une raideur diff£rente de celle correspondant k 
la loi de comport ement . 

Selon une variante, l'etape d' enregistrement du 
modele numerique de 1' ensemble d' articulations standard 

20 consiste a definir une alternance de corps rigides et 
d' articulations , et a definir pour chacun des couples de 
corps un ensemble de parametres num^riques caracterisant la 
mobilite ou raideur globale resultant de 1 7 action de 
1' ensemble des elements intercalaires [par exemple disques 

25 intervertebraux] et des elements de liaison [par exemple 
ligaments] ayant un effet sur les parametres d' interaction 
[rigidites] entre les deux corps. 

Selon un mode de mise en ceuvre particulier, l'6tape 
de personnalisation consiste a acquerir au moins une image 

30 de 1' ensemble d' articulations d'un patient donne, a 
extraire de ladite image les informations necessaires a la 
construction d'un modele reel par reconnaissance de la 
position des articulations visibles dans ladite image, et 
de modifier le moddle standard en fonction dudit module 

35 reel . 
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Avantageusement, I'etape d'enregistrement d'un 
modele numerique consiste a definir un ensemble standard de 
donnSes numeriques comprenant pour chacune des 
articulations representees sous la forme d'un corps 
5 rigide : 

• un premier descripteur giometrique de 
position de reference, correspondant a la geom^trie 
de 1' ensemble d' art iculat ion pour un patient 
« standard » dans une position «de ref6rence», 

10 ledit descripteur §tant determine pour chaque corps 

rigide de fagon relative par rapport a un corps 
adjacent ; 

• un deuxi erne descripteur mecanique 
d' interaction avec chacun des corps adjacents, 

15 ledit descripteur mecanique etant representatif de 

la loi de comportement lorsque 1' ensemble 
d # articulation est soumis a I 7 action d'au moins une 
contrainte exterieure ; 

20 I'etape de personnalisation consistant a modifier 

ledit ensemble standard de donnees par des donnees 
personnalisees . 

Selon un mode de mise en oeuvre particulier, le 
25 procede comporte en outre une 6tape de recalage consistant 
a met t re en correspondance des donnees d' images radios et 
des donnees d' acquisition externe, cette etape se 
decomposant en deux sous -Stapes : 

• recalage de la reconstruction radio par 
30 rapport a la courbe 3D issue de donnees 

d' acquisition externe dans la mime 
position ; 

• determination de la repartition des points 
de la courbe 3D associes aux vertdbres, 



WO 2004/061721 



9 



PCT7FR2003/003943 



posit ionnes dans le repere de Stokes et leur 
tangente associee. 



On comprendra mieux 1' invention a l'aide de la 
5 description, faite ci-apres a titre purement explicatif, 
d'un mode de realisation de 1' invention, en reference aux 
figures annexe es : 

• la figure 1 illustre 1 ' architecture de la 
personnalisation mecanique du modele ; 

10 • la figure 2 represente 1 ' architecture generale du 

simulateur selon un mode de realisation de 
1' invention ; 

• la figure 3 represente 1 ' architecture du modele 
selon un mode de mise en ceuvre de 1' invention ; 

15 • la figure 4 illustre la mise en oeuvre du proc^de 

de recalage ; 

• la figure 5 illustre la definition des angles 
intervertebraux ; 

• la figure 6 represente un exemple des relations 
20 entre les angles intervertebraux dans le plan 

frontal (abscisse en ° et ordonn^e en °) ; 

• la figure 7 illustre le calcul des centres de 
rotation ; 

• la figure 8 represente un exemple de calcul des 
25 centres de rotation sur un sujet scoliotique ; 

• la figure 9 illustre 1 ' interpolation des 
couronnes ; et 

• la figure 10 represente un module de repartition 
des masses dans le tronc. 

30 

L' invention se distingue des produits connus de 
guidage chirurgicaux, utilises pour assister le chirurgien 
pendant son operation. En effet, 1' invention consiste 
notamment a modeiiser la colonne vertebrale de l'individu 
35 qui va etre oper£, de simuler la pose de 1' implant ou de 
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prothdse et de calculer la position d'^quilibre de 
1'individu une fois 1' implant ou la prothdse pose . 

Le traitement des donnees etant logiciel, 
1' invention concerne aussi 1 ' architecture logicielle mise 
5 en place pour la realisation de 1 ' architecture 
fonctionnelle. La structure logicielle comprend plusieurs 
serveurs de bases de donnees : une base de donn§es 
utilisateurs, une base de donnees vert^bres, une base de 
donn§es patients et une base de donnees implants . Les 
10 donnees contenues dans ces bases sont interrog^es et mises 
a jour par diff^rents utilisateurs afin de construire le 
modele 3D de la colonne du patient puis de simuler les 
consequences de la mise en place d'un implant. 

Le systdme selon 1' invention est un simulateur des 
15 consequences biomecaniques et cin^matiques de l'acte 
chirurgical des traitements des pathologies rachidiennes 
sur la morphologie du patient. La Figure 2 represente 
1 ' architecture generale du simulateur. 

II doit permettre & un chirurgien d'optimiser et 
20 ameliorer la planif ication de son traitement. Ce simulateur 
pourra lui apporter une meilleure connaissance des 
propri^tes g^ometriques et mecaniques des divers tissus de 
la colonne vertebrale. II pourra aussi tester differentes 
approches de son geste pour permettre une correction 
25 optimale . 

II doit r^pondre & un besoin de sante, puisque la 
tendance est k la recherche du confort, de la security, de 
la qualite et de la fiabilit£. L'alldgement des traitements 
m^dicaux se traduit par une meilleure quality de vie post- 
30 operatoire et en particulier par des hospitalisations de 
frequence moindre et de duree plus courte . 

Enfin, il peut egalement avoir des retombees dans 
le domaine de 1 1 enseignement pour 1 ' apprentissage des 
futurs chirurgiens. 
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Le simulateur est un outil d'aide a 1' analyse de la 
faisabilite du geste. II permet de simuler l'6quilibre pre- 
operatoire, les efforts inter-segmentaires correspondants 
tenant compte de l'effet des rigidit^s musculaires et 

5 ligamentaires, apprehender 1' evolution de cet equilibre et 
de ces efforts en post-operatoire en fonction des courbures 
introduites par le chirurgien lors de 1' operation. Le 
chirurgien devra indiquer en entree la deformee voulue de 
la zone instrument§e. Pour cela une approche dite 

10 « globale » sera utilisee. 

Deux radiographies sont necessaires pour la 
realisation d'une simulation. Les images £ ajouter au 
systeme peuvent §tre : 

- le fichier r§sultat de la numerisation (scan) 
15 d'une radiographic traditionnelle ; 

- un fichier fourni par un autre appareil de 
radiologie (radiographic numerique) . 

L' alimentation du systeme peut §tre manuelle (un 
utilisateur va fournir les fichiers d' images au systeme) ou 
20 automat ique, les images sont stock^es directement par les 
appareils de radiologie et r§cuperees via un reseau 
intranet/ internet . 



Deux sortes de numerisations seront possibles : 
25 - numerisation manuelle : des points specifiques 

doivent etre reperes manuellement par un acteur sur ces 
images ; 

- numerisation semi -automat ique . 



30 Les points specifiques detectes par la numerisation 

vont permettre de calculer les coordonn^es 3D des 
vertebres (donnees geom^triques) . 

Les donnees geom^triques issues des images 
num^riques vont permettre de construire un modele de la 

35 colonne du patient en trois dimensions. Ce modele resulte 
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de 1' adaptation d'un module 3D standard, pr^defini dans le 
syst^me, aux caracteristiques g6om£triques du patient. 

L' utilisateur doit pouvoir visualiser le modele 3D 
du rachis selon les plans frontal, sagittal et apical ; de 
5 plus, il doit pouvoir le comparer aux radiographies 
utilis€es pour sa construction. 

Les caracteristiques mecaniques du patient 
(rlsultats issus d' acquisitions (tests cliniques) ) 
permettent de personnal iser le modele. Pour cette 
10 operation, un modele g^ometrique standard de la colonne 
vertebrale est enrichi par les donn6es mecaniques 
personnelles du patient concerne. 

Les caract6ristiques biomecaniques de la colonne du 
patient (angle de scoliose, rotation axiale, etc.) sont 
15 calculees a partir des radiographies du patient (avant la 
premiere simulation) , ou a partir du module resultant d'une 
simulation. Certains parametres pourront etre pr^-calcules 
(pente sacree, angle d' incidence, gite sagittale, courbures 
rachidiennes) sur les radiographies. 
20 Lors de la simulation, 1 ' utilisateur doit avoir la 

possibility de choisir des segments de colonne sur lesquels 
il impose des deplacements . Lors de la simulation de 
l'6quilibre, 1 ' utilisateur va pouvoir entrer les valeurs de 
certains parametres cliniques qu'il souhaite simuler pour 
25 un segment donn£ . L' utilisateur doit pouvoir consulter les 
mobilites ou rigidites du rachis, la graduation sera 
normee. Le type de la tige utilisee pour deformer un 
segment peut etre choisi parmi plusieurs propositions 
(rigidity, diamdtre, etc.). 
30 A la suite d'une manipulation du modele par 

1' utilisateur, le logiciel doit verifier la validite des 
actions effectuees, et prevenir 1 ' utilisateur de toute 
incoherence (valeur impossible a obtenir) . 

Le simulateur doit afficher les nouvelles courbures 
35 du modele et la nouvelle position (stature) de l'gquilibre 



WO 2004/061721 



13 



PCT/FR2003/003943 



du patient. II doit etre possible de comparer la (les) 
courbure(s) avec la (les) courbure(s) initiale(s) (avec le 
modele 3D ou la radiographie) . 

L' evolution relative des efforts intervertebraux 
5 est quantifiee et representee graph i quement . 

L' utilisateur doit pouvoir comparer les rigidites 
du rachis avant et aprds la simulation. 

Les efforts resultants dans la tige suite aux 
manipulations effectuees doivent indiguer si la tige va se 
10 deformer. 

Des dossiers patients, contenant les informations 
concernant la simulation, doivent pouvoir §tre crees par un 
utilisateur et maintenus a jour par le systeme, (sauvegarde 
de plusieurs simulations dans un meme dossier) . De plus, le 
15 systeme doit pouvoir retirer ou saisir les informations 
concernant un patient (Sge, poids, taille, etc.) dans des 
systemes de donnees externes, sous reserve de leur 
existence . 

Les droits et profils utilisateurs vont permettre 
20 de differencier les f onctionnalites disponibles pour chaque 
utilisateur. En fonction de leurs droits et profils, les 
utilisateurs auront acces a des interfaces graphiques 
differentes. De plus, un systeme de preference pourra etre 
mis en place pour chaque utilisateur (position des menus, 
25 page d'ouverture etc.), ainsi qu'un mecanisme de 
tragabilite, qui permettra de suivre les evolutions d'un 
dossier patient. 

A tout moment, et en fonction de ses droits, un 
utilisateur peut avoir acces aux donnees d' entrees, c'est- 
30 a-dire visualiser, remplacer ou modifier les images 
utilisees pour la construction du modele, les points 
numerises ou les operations effectuees pour le traitement 
semi- automat ique de ces images, les implants places sur le 
module, les caracteristiques mecaniques du patient... 
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L'historique d'une simulation represente la 
totalite des actions effectives sur un patient (parametres 
de la simulation) . L' utilisateur doit pouvoir revenir sur 
une action passee (changer une radiographic, ressaisir les 
5 points pour la numerisation, calculer de nouveaux 
parametres cliniques, ou modifier les actions faites au 
cours de la simulation) 

La Figure 1 represente 1 ' architecture de la 
personnalisation mecanique du modele. 
10 La personnalisation mecanique du modele est basee 

sur trois types de donnees : 

- les radiographics du patient avec des reperes 
cutanea ; 

- 1' acquisition de la courbure generale du rachis 
15 sous diff^rentes postures caracteristiques ; 

- des donnees anthropometriques . 



A partir de ces donnees # des traitements utilisant 
les lois de la mecanique nous permettent d'obtenir : 
20 - la g6om£trie du patient ; 

- les parametres cliniques ; 

- le modele mecanique personnalis£ . 

La numerisation des radiographics du patient nous 
25 permet d'obtenir une precision de six points par vertebre . 
En effectuant une extrapolation, on arrive a une precision 
de douze points par vertebre. 



L'^tape suivante est 1 ' acquisition de la courbure 
30 generale du rachis sous diff^rentes postures 
caracteristiques . 

Pour ce faire, lors d'un examen clinique, le 
patient etudie est soumis a une s^rie de tests au cours 
desquels sera §valu£e la ligne generale de la colonne au 
35 niveau dorsal. Le positionnement du bassin et des epaules 
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sera necessaire pour d^finir les orientations des 
extremit^s de la colonne. 

Le patient devra §tre, eventuellement, maintenu au 
niveau du bassin pour limiter 1 ' intervention des membres 

5 exterieurs dans 1 ' etablissement de son £quilibre general, 
il sera effectue une acquisition i. 1'aide d'un appareil, 
des positions dans 1'espace de reperes cutanes 
identif iables lies aux vertibres scoliotiques par rapport a 
un referentiel connu. 

10 Une acquisition sera effectuee patient au repos 

dans des conditions aussi proches que possible des 
conditions de prise de la radiographic calibr^e patient 
debout au repos. Cette acquisition sera utilisee pour 
determiner la ligne des centres des corps vertebraux £ 

15 partir de la ligne generale de la colonne. Ceci est rendu 
possible par le fait que la radio calibree peut dormer les 
coordonn§es de tous ces points grace aux reperes plombes. 
Le positionnement des vertebres peut done etre determine en 
fonction des reperes cutanes pour une position donnie , 

20 debout, au repos. II sera ensuite etabli une transformation 
corrigee entre les reperes cutanes et le positionnement des 
vertebres prenant en compte 1' influence de la cin^matique 
des vertebres entre elles lors des mouvements du patient. 

25 A partir des donnles anthropometriques, on 

determine les masses segmentaires et les centres de masses, 
permettant ensuite de determiner les moments et les centres 
de masse par niveau vertebral. 

30 L ' action musculaire sera dissociee de 1' action 

inter segmentaire . 



35 



Les donn6es disponibles dans la litterature 
permettent d'evaluer qualitativement la forme des lois de 
comportement admissibles pour le moddle. 



WO 2004/061721 



16 



PCT/FR2003/003943 



Les lois de comportement doivent repondre aux 
exigences suivantes : 
5 la plupart des lois doivent suivre des 

comportements impairs ; 

- le comportement asymptotique doit etre assure ; 

- les phenomenes de couplage sont pris en comptes ; 

- la pertinence et la simplicity de calcul . 

10 

Les lois de comportement devront etre recalculees 
pour chaque niveau vertebral et pour chaque patient af in de 
prendre en compte les singularites dues aux pathologies 
etudiees . 

15 

Selon une variante de 1' invention, on met en 
correspondance les radios avec les acquisitions externes. 
La mise en correspondance se passe en deux Stapes : 

• Le recalage de la reconstruction radio par rapport k 
20 la courbe 3D issue de donnes d' acquisition externe 

dans la meme position. 

• La determination de la repartition des points de la 
courbe 3D associes aux vertebres, positionn^s dans le 
repere de Stockes et leur tangente associees. 

25 Lors des acquisitions radio et de donnees externes, trois 
reperes spatiaux sont §galement enregistres. Les deux 
extr6mites de la rSglette, et une bille en plomb (point 
bas) , placee en bas de la colonne. 

30 Ces trois points dont on connait la position 

relative, a la fois dans le repere de reconstruction radio 
et dans le repere des donnees d' acquisition externes, 
permettent de recaler les deux enregistrements dans le meme 
repere, celui de la radio. 
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Les trois points sont suffisants a determiner 
un repere commun dans lequel 1' ensemble des donnSes est 
connu, ce qui permet done de definir la matrice de passage 
entre les deux enregistrements . 
5 Compte tenu des erreurs d' acquisition possible, 

rotation du poignet involontaire, rugosite de la peau 
variable le long de 1 ' enregistrement , trois autres 
recalages sont egalement effectues : 

1) on compare la distance entre le point 
10 bas de la reglette et le point bas plac6 sur la peau, 

pour la radio et 1 ' acquisition externe . La difference 
entre les deux enregistrements quantifie le trajet 
effectue par la spline sur la peau et done reajuste 
1 ' acquisition vis-a-vis de la pression exercee par 
15 1'operateur sur la peau du sujet. 

2) Cet ajustement est complete par un 
ajustement des longueurs des acquisitions externes de 
telle mani^re a garantir le fait que les quatre 
acquisitions partent du m§me point. 

20 3) Lors de 1 ' acquisition, suivant le 

positionnement de 1'operateur, il est frequent de 
constater un biaisement de 1 ' acquisition vers la 
gauche ou la droite. Ce decalage est corrig6 en 
redressant les acquisitions en contraignant 

25 1' acquisition allant jusqu'au point bas de respecter 

la jonction entre le point concourant d^fini au 2 et 
le point bas d£fini sur la radio. 

Suite a ces trois recalages, il est possible 
d' avoir une confiance certaine dans le recalage effectue, 

30 corrig6 des erreurs. 



Repartition des vertebres sur la courbe 3D : 
A partir du recalage entre les radios, la 
deuxieme #tape consiste a placer les vertebres sur la 
35 spline. On considere que le point de correspondance entre 
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la vertebre et la spline est le point d' intersection du 
repere de Stockes (plan XY) avec celle-ci. 



Spur les vertlbres du bas de la colonne, la 
5 distance importante entre la colonne et la surface du dos, 
ainsi que la forte inclinaison des vertebres provoquent des 
positionnements non successifs des vertebres sur la Courbe 
3D. Sur cette portion d 1 enregistrement, on applique la 
demarche inverse, on place des points 6quidistants puis on 
10 definit leurs positions dans le repdre de Stokes de la 
vertebre correspondante . 

On enregistre ggalement les tangentes de la 
Courbe 3D aux points consideres de maniere a pouvoir par ia 
suite repositionner le modele. 
15 Cette methode permet de d^finir egalement sur 

les enregistrements d' acquisitions externes la position du 
sacrum et 1' angle du bassin avec la verticale, appele angle 
sacre . 

Meme si ces parties anatomiques ne font pas 
20 parties de la reconstruction k partir de la radio, ils sont . 
utiles pour definir l'€quilibre du rachis. 

On suppose ensuite que la repartition des 
vertebres sur toutes les acquisitions externes reste la 
m g me sur toutes les acquisitions, ce qui permet de 
25 positionner les vertebres dans toutes les positions. 



Notons quelques remarques sur 1' analyse des donn6es : 

o L ' etude se fait sur les deux plans 
30 frontal et sagittal separement . En effet, le choix a 

ete pris d'etudier dans un premier temps la colonne 
vertebrale de manidre d^couplee. C'est pourquoi les 
mouvements demandes aux patients sont des inflexions 
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laterales pour 1' analyse du plan frontal et des 
flexions/extensions pour 1' analyse du plan sagittal. 

Le calcul des angles intervertebraux est 
5 illustre sur la figure 5. 

Le but est de trouver une relation liant les 
angles intervertebraux. Pour cela, un module d' analyse a 
ete ajoute, tragant un angle intervertebral en fonction de 

10 celui qui se trouve « en dessous » lorsque l'on parcourt la 
colonne de bas en haut . 

Exemple de trac6 des relations entre les angles 
intervertebraux (le premier graphique reprlsente 0(C7,n)/<n,r2) 
haut a gauche, puis on parcourt la colonne en descendant 

15 dans le sens de lecture classique) :La figure 6 illustre un 
exemple des relations entre les angles intervertebraux dans 
le plan frontal (abscisse en ° ; ordonnee en °) 

Les conclusions de ces nuages de points sont 
20 les memes pour tous les sujets. La distribution est 

consider6e lineaire. Le module est done batit sur ces 

constatations : 

On definit 1' equation de regression sur les 

points a chaque niveau intervertebral . A chaque equation on 
25 adjoint des bornes caract6risant la nature born£e de la 

mobilite d'une vertebre par rapport a celle qui se trouve 

en dessous. On fait alors l'hypothese que les mouvements 

d' inflexion lat^rale et de flexion/extensions maximales 

d£finissent ces but^es. 

30 

Cette approche a £te validee en analysant un 
certain nombre de sujets sains. Pour la plupart des sujets, 
et des niveaux intervertebraux, la mod^lisation adoptee 
semble tout a fait appropriee : les coefficients de 
35 correlation sont le plus souvent super ieurs a 0,8. 
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Cependant , il existe des niveaux pour lesquels la 
modelisation est incorrecte. Ce ph^nomene correspond k des 
niveaux bloques ou presque bloques chez le patient. On 
considere done qu'il reste correct de construire la 
5 regression lineaire, et les but£es contraindront de ce fait 
suf f isamment le modele. 

Cette methode permet de definir des lois de 
comportement du disque T1/T2 k L5/S1. 

10 Les centres de rotations d'une vertebre par 

rapport k une autre se calculent a partir des differentes 
positions enregistr^es dans le plan frontal (bending) et 
sagittal (Flexion/extension) . Les points de la vertebre 
supirieure correspondant a chaque position sont places dans 

15 le repere de Stockes de la vertebre inferieure. Ces points 
forment alors une trajectoire consid6ree comme circulaire 
dont on calcule le centre de rotation par la methode des 
moindres carres . 

Ceci est illustre sur la figure 7. 

20 

Cette methode permet d' identifier des centres 
de rotation apparents et personnal is§s verifiant les 
donnees experimentales obtenues sur un patient. 

Un exemple de calcul des centres de rotation 
25 sur un sujet scoliotique est illustre sur la figure 8. 

On constate que cette methode permet de 
retrouver la lordose lombaire a partir des acquisitions 
externes . Cette methode necessite un minimum de trois 
30 acquisitions distinctes. 

Ces centres de rotations sont calcules dans le 
repere global de la radio. Pour chaque vertebre, deux 
centres de rotation distincts sont calcules dans l'espace, 
un pour les mouvements sagittaux et un pour les mouvements 
35 frontaux. Ces centres sont ensuite exprimes dans le repere 
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de Stokes de la vertebre inferieure de maniere a assurer le 
repositionnement dans le solver. 



Un exemple d' algorithme est repr^senti sur la 
5 page suivante. 
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• Acquisition du modfele cindmatique : 

Module geometrique issue des radios 

Relation Unmake entre les differents angles intervertebraux 

Amplitude maximum des differents angles intervertebraux 

Position dans les repfcres de Stokes des vert&bres des centres de rotation frontaux et sagitta 



• Initialisation du variant : angle intervertebral L5/S 1 

• Calcul des angles intervertebraux initiaux sur le modele geometrique 

Tant que le critfere de confort n'est pas optimise 
Faire 

L'angle variant est incr&nentd du pas de calcul (par defaut 0,001) 

l Si Tangle intervertebral variant ddpasse ses butees Alors Tangle vaut la valeur de la butee 
et le variant est Tangle suivant 

| Pour tous les angles suivants 
Faire 

L' angle intervertebral i est incr^mente en fonction de Tangle i+1 

Si Tangle intervertebral i depasse ses butfes Alors Tangle vaut la valeur de la butee 
Fin Faire 

Repositionnement du module k partir des nouveaux angles intervertebraux. 
I Application du critere mecanique pour repositionnement du module autour de T axe des 
tStes fSmorales 

Recalcule des angles intervertebraux sur le modele geometrique suite & sa deformation 
Calcul du critere de confort 

Si le niveau de confort diminue, on incremente le modele dans T autre sens. 

Si les deux sens de parcours sont faits, on diminue le niveau d* exigence du critere de 

confort 



Fin Faire 
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La colonne vertebrale est jusqu'3. maintenant 
simulee cin^matiquement. On aboutit a une geometrie du 
rachis repr^sentant l'equilibre du patient. Il est 
5 interessant dgsormais de connaitre la repartition des 
efforts qui s'exercent sur la colonne, afin de construire 
un modele mecanique complet. L'art anterieur connait une 
approche pour le calcul des efforts s'inspirant du modele 
anthropom^trique . Le principe de cette approche est de 

10 decouper le tronc en quatre tranches rattachees a quatre 
parties bien precises de la colonne vertebrale et de 
calculer leur poids et centre de gravite. Le contour de 
chaque tranche est obtenu par dimensionnement de patron 
g^nerique de la population moyenne aux dimensions du sujet 

15 mesure . 

La m€thode utilis^e dans le module selon 
1' invention est sensiblement differente. Cinq couronnes 
sont enregistrees & partir d' enregistrements de contours. 
Ces couronnes sont ensuite discr£tis£es en 60 points. Dans 

20 un deuxiSme temps, ces couronnes sont positionnees par 
rapport & la colonne §l partir du recalage. Ce 
positionnement se fait en confondant le point sur la peau, 
au niveau de la vertebre associee & la couronne, enregistre 
par la courbe de recalage avec le point de la couronne 

25 correspondant au centre du dos . Les autres couronnes sont 
ensuite interpolees pour obtenir 18 couronnes definissant 
17 tranches centrees verticalement sur le centre du corps 
vertebral de la vertebre correspondante . 

30 

La Figure 9 illustre 1 ' interpolation des couronnes. 
Le volume et le centre de gravite de chaque tranche ainsi 
definie est calcule en les decoupant en prismes 
61ementaires . 
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A partir des volumes et des centres de gravit^s de chaque 
tranche, on dresse le modele suivant : 



5 ■=> Le poids de la tSte et des bras est repris uniquement 

par la colonne vertebrale 
«=> Chaque tranche est divis^e en deux : les visceres, et 
les parties solides comme les cotes, muscles, peau 

10 o La fraction solide de la tranche est reprise 

entierement par la colonne 
o La fraction viscere a un comport ement 
hydrostatique : 

■ Les visceres n'engendrent pas de moment 

15 ■ Seule la composante normale a la colonne 

vertebrale voit les efforts transmis par 
les visceres 

■ Les composantes tangentielles a la colonne 
sont transmises a la tranche d'au-dessous 

20 

La Figure 10 repr^sente un modele de 
repartition des masses dans le tronc . 

On obtient alors le torseur des efforts sur chaque vertebre 
dependant uniquement de la geom£trie entree . 

25 La repartition viscere / corps dur pour chaque tranche est 
pour 1' instant d£fini & partir d'un nombre restreint de 
donnees . Cette repartition sera par la suite affinee. 

L' architecture du programme se base sur une 
structure modulaire correspondant a chaque etape de calcul 

30 independant, de maniere a pouvoir faire evoluer le 
programme sans alterer son f onctionnement . Les modules de 
calcul sont presentes ci-dessous. lis sont classes en 4 
categories : 

• Les modules d' acquisitions 

35 • Les modules de definition du modele cinematique 



WO 2004/061721 



25 



PCT/FR2003/003943 



• Les modules du solveur de l'equilibre 

• Module de calcul des efforts sur le rachis 
instruments. 

L' invention est decrite dans ce qui precede a titre 
d'exemple. II est entendu que l'homme du metier est £ meme 
de realiser differentes variantes de 1 ' invention sans pour 
autant sortir du cadre du brevet. 
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REVINDICATIONS 



1 - Procede de simulation biom^canique d'un 
ensemble d' articulations osseuses d'un patient, notamment 
5 du rachis, comportant : 

• une etape d'enregistrement d'un modele 
numerique tridimensional representee au moins en 
partie par des corps rigides relies par des 
articulations, dans une position de reference ; 

10 • une etape de personnalisation de la 

g^ometrie dudit modele [position relative dans 
l'espace de chacun desdits corps rigides] par des 
donnees specifiques k un patient dans ladite 
position de reference [par exemple des 

15 radiographics du patient] ; 

• une etape de personnalisation dudit modele 
numerique par part icularisat ion des parametres 
d' interaction [de mobilites ou caracteristiques de 
rigidite] de chacune des articulations reliant 

20 lesdits corps rigides en fonction des 

caracteristiques constatees sur le patient ; 



caracterise en ce que 



25 1' etape de part icularisat ion des parametres 

d ; interaction [de mobilites ou caracteristiques de 
rigidite] consiste a : 



acquerir les positions dans l'espace d'une 
30 partie au moins des corps rigides, et a proceder a 

une interpolation pour determiner la position 
calcuiee des autres corps rigides pour construire 
une table numerique comportant les positions 
relatives de chacun des corps rigides ; 
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• exercer au moins une contrainte determinee 
sur le patient et £ acqu^rir une information sur la 
position generale d'equilibre resultante du patient 
[autre que la position de reference] ; 
5 • determiner des fonctions analytiques 

permettant d'approximer les paramdtres 
d' interaction [de mobilites ou caract^ristiques de 
rigidite] afin de reproduire les positions 
relatives mesur6es [fonctions desdites contraintes 
10 appliqu^es sur le patient et des modifications 

geom^triques observes sur le patient, et resultant 
de ces contraintes] et ce pour chaque couple de 
corps rigide. 

15 

2 - Procede de simulation biomecanique d'un 
ensemble d' articulations osseuses selon la revendication 1, 
caracterise en ce que le module numerique est defini par 
des parametres de positions g^ometriques des corps rigides 

20 et par des parametres de rigidity des articulations reliant 
les corps rigides . 

3 - Procede de simulation biomecanique d'un 
ensemble d' articulations osseuses selon la revendication 1 

25 ou 2, caracterise en ce que l'etape de representation du 
resultat d'une contrainte consiste a recalculer le module 
personnalis§ [en position d'equilibre] resultant d'un 
ensemble de contraintes [par exemple implantation d'une 
prothese ou d'un implant] comprenant au moins une 

30 contrainte statique exercle sur au moins deux corps 
rigides, et imposant un positionnement relatif avec une 
mobilite ou une raideur differente de celle correspondant a 
la loi de comportement . 
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4 - Proced§ de simulation biomecanique d'un 
ensemble d' articulations osseuses selon les revendications 
1 et 2, caract6ris6 en ce que l'etape d'enregistrement du 
moddle nurmirique de 1' ensemble d' articulations standard 

5 consiste a definir une alternance de corps rigides et 
d' articulations, et a definir pour chacun des couples de 
corps un ensemble de paramltres num#riques caracterisant la 
. mobilite ou raideur globale resultant de 1' action de 
1' ensemble des elements intercalaires [par exemple disques 
10 intervertebraux] et des Elements de liaison [par exemple 
ligaments] ayant un ef f et sur les paramdtres d' interaction 
[rigidit^s] entre les deux corps. 

5 - Precede de simulation biomecanique d'un 
15 ensemble d' articulations osseuses selon 1'une au moins des 

revendications pr<§c§dentes , caracteris£ en ce que l'etape 
de personnalisation consiste £ acquerir au moins une image 
de 1' ensemble d' art iculations d'un patient donn£, a 
extraire de ladite image les informations necessaires £ la 
20 construction d'un modele reel par reconnaissance de la 
position des articulations visibles dans ladite image, et 
de modifier le moddle standard en fonction dudit modele 
reel. 

25 6 Procede de simulation biomecanique d'un 

ensemble d' articulations osseuses selon l'une au moins des 
revendications precedentes, caract£rise en ce que l'etape 
d' enregistrement d'un modele numerique consiste a definir 
" un ensemble standard de donnees numeriques comprenant pour 
30 chacune des articulations representees sous la forme d'un 
corps rigide : 

• un premier descripteur geometrique de 
position de reference, correspondant £ la g^ometrie 
de 1' ensemble d' art iculat ion pour un patient 
35 « standard » dans une position «de references, 
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ledit descripteur etant determine pour chaque corps 
rigide de fagon relative par rapport a un corps 
ad j acent ; 

• un deuxieme descripteur mlcanique 
5 d' interaction avec chacun des corps adjacents, 

ledit descripteur mecanique £tant representatif de 
la loi de comportement lorsque 1' ensemble 
d' articulation est soumis a 1' action d'au moins une 
contrainte exterieure ; 
10 1' etape de personnalisat ion consistant a modifier ledit 

ensemble standard de donnees par des donnees 

personnalisees . 

7 - Proc§d6 de simulation biomecanique d'un 
15 ensemble d' articulations osseuses selon l'une au moins des 
revendicat ions precedentes, caracterise en ce qu'il 
comporte en outre une etape de recalage consistant a mettre 
en correspondance des donnees d' images radios et des 
donnees d' acquisition externe, cette etape se decomposant 
20 en deux sous -Stapes : 

• recalage de la reconstruction radio par 
rapport a la courbe 3D issue de donnees 
d' acquisition externe dans la meme 
position ; 

25 • determination de la repartition des points 

de la courbe 3D associes aux vertebres 7 
positionnes dans le repere de Stokes et leur 
tangente associee. 
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Figure 1 

Architecture de la personnalisation mecanique du modele 
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Figure 2 
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Figure 5 
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Figure 6 
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Figure 7 
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Figure 8 
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Figure 9 



•1*? 



a= 



i 




10/141250 

WO 2004/061721 PCT/FR2003/003943 

10/10 

Figure 10 




